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Abstract 

Measuring external sources of background for a deep underground laboratory at the 
Homestake Mine is an important step for the planned low-background experiments. 
The naturally occurring 7-ray fluxes at different levels in the Homestake Mine are 
studied using Nal detectors and Monte Carlo simulations. A simple algorithm is 
developed to convert the measured 7-ray rates into 7-ray fluxes. A good agreement 
between the measured and simulated 7-ray fluxes is achieved with the knowledge of 
the chemical composition and radioactivity levels in the rock. The neutron fluxes 
and 7-ray fluxes are predicted by Monte Carlo simulations for different levels includ- 
ing inaccessible levels that are under construction for the planned low background 
experiments. 

Key words: 7-ray flux, neutron flux, Nal detectors, Monte Carlo simulation, 
underground laboratory, radioactive decays, dark matter, double-beta decay 
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1 Introduction 



Neutrinos and dark matter are believed to hold the key to physics beyond the 
Standard Model [T|2|3|^f5|6|7] . However, very little is known about the gen- 
eral properties of neutrinos such as their absolute mass and magnetic moment, 
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whether they are Dirac or Majorana particles, or the number of species. Even 
less is known about the nature and quantity of dark matter in the universe. 
A deeper understanding of both neutrinos and dark matter is important to 
understand physics beyond the Standard Model. Planned underground low- 
background experiments are intended to probe these properties using rare 
event physics processes such as neutrinoless double-beta decay, neutrino oscil- 
lations, and the direct detection of dark matter. Because of a very low event 
rate, the reduction of backgrounds from radioactivity and cosmic rays is es- 
sential to these experiments. 

The solution to reducing cosmic rays is to build experiments in a deep under- 
ground laboratory. However, deep underground experiments require appropri- 
ate shielding against 7 rays and neutrons from natural radioactivities in the 
surrounding materials and rocks. As the major contributors to the radiation, 
the 7 rays and neutrons range from keVs to MeVs in kinetic energy. The 7 
rays are mainly from the decays of 238 U, 232 Th, and 40 K in the rocks. The 
neutrons come primarily from (a,n) reactions induced by 238 U, and 232 Th de- 
cays and fission decays of 238 U [8f9] . Care must be taken to understand the 
7-ray and neutron fluxes as a function of energy in order to design appropriate 
shielding. This is usually done by measurements in combination with Monte 
Carlo simulations. SNO [ID], Gran Sasso [TT|T3j .Frejus [H], Boulby [15], and 
other underground laboratories [16] have reported respective 7-ray and neu- 
tron fluxes in terms of the measurements and Monte Carlo simulations. 

The Homestake Mine in Lead, South Dakota, has been selected by the Na- 
tional Science Foundation (NSF) as the potential site for Deep Underground 
Science and Engineering Laboratory (DUSEL). As funding for DUSEL is being 
secured, an early science program is being developed at the interim Sanford 
Laboratory. This makes the characterization of the 7 rays and neutrons in the 
Homestake Mine a necessary endeavor. Generally, the fluxes of 7 ray and neu- 
trons reveal the concentrations of the radioactive sources and the composition 
of the surrounding rocks. Because the geology of Homestake Mine varies from 
location to location, the variations in the radioactivity concentration in the 
rock [21] result in site-dependent 7-ray and neutron fluxes. 

We demonstrate for the first time the site-dependent 7-ray flux with measure- 
ments at different levels using Nal detectors in Section 2. We measure the 
7-ray fluxes at the surface in the administration building, the underground 
800-ft level, 2000-ft level, and 4550-ft level. Several measurements at differ- 
ent spots separated by no more than 100 feet were done at the same level. 
Since many rock samples in the different levels at Homestake Mine have been 
collected and analyzed [20], we are able to perform a Monte Carlo simula- 
tion to estimate the 7-ray fluxes for different levels. The simulated results are 
discussed and compared to the measurements in Section 3. 
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The neutron fluxes at different levels are characterized through a Monte Carlo 
simulation. The neutron yield and energy spectrum in the rock are calcu- 
lated utilizing the tools and methods described in Ref. [8]. Those neutrons 
are tracked through the rock to the experimental hall. We elaborate on the 
neutron fluxes in Section 4. 



2 The 7-ray fluxes at different levels 



The 7 rays from natural radioactivity in the rocks at the Homestake Mine are 
measured using Nal detectors. Three identical Nal detectors were deployed 
to directly measure 7 rays at different levels. The Nal detectors are 3 x 3 
from Bicron Crystal Corporation. All three detectors are read out using a 
Digital Gamma Finder Four Channel (DGF4C) CAMAC module supplied by 
X-Ray Instrumentation Associates (XIA) [17] . This module is a 14-bit digitizer 
with 40 MHz sampling rate. The CAMAC create is connected to an Express 
card of a Hewlett Packard (model 8730) laptop running Microsoft Windows 
Vista. The system was controlled using the standard software provided by 
XIA. This data acquisition software runs in the IGOR Pro environment [TS] 
and produces binary data files that are read in and analyzed using the ROOT 
framework [T9], 

Three detectors were arranged horizontally in the same plane and were sepa- 
rated by about 50 cm from each other. The measured 7-rays energy spectrum 
from each detector provides a cross-check to spot anomalous detector response, 
which can then be corrected on the site. The internal background from the 
detectors was measured at the surface lab with lead bricks that surround the 
detectors. The contribution from the internal background to the energy range 
of interest (less than 3 MeV) is less than 0.1% of the measured external 7 
rays. 

The detectors were calibrated using a 60 Co 7-ray source. This simple calibra- 
tion allows us to understand the energy response and identify the well-known 
peaks from the external radiation. The four well-known peaks, 0.609 MeV and 
1.764 MeV from 214 Bi in 238 U decay chain, 1.461 MeV from 40 K decay, and 
2.615 MeV from 232 Th decay chain, are distinguished. These four peaks are 
used as the calibration points for each collected data file. Despite the small 
variations in the gain of the detectors resulting in an energy shift in the spec- 
trum, the four well-known peaks always provide a reliable self-calibration to 
the entire measured spectrum. 
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y-Ray Energy (keV) 



Fig. 1. Shown is the measured differential 7-ray flux at the surface administration 
building and underground on the 800-ft, 2000-ft, and 4550-ft levels from one of the 
three detectors. 

The differential flux is determined by using the formula 

where $(£') is in unit of cm _2 s _1 keV _1 , R is the registered count rate (counts/s) 
of the detector for a selected energy region, cr(E) is the photon attenuation 
factor (cm 2 /g) of Nal with respect to the selected energy [22], and m is the 
effective mass (g) of the Nal crystal. The term cr(E) (cm 2 /g) • m(g) accounts 
for the effective area of the detector that subtends the 7 rays and the detec- 
tion efficiency of the 7 rays as a function of energy. The value of this term is 
obtained by the combination of the calibration and Monte Carlo simulation. 

The final experimental results for the surface administration building and 
underground on the 800-ft, 2000-ft and 4550-ft levels are shown in FigureU 
The 7-ray flux with energy greater than 0.1 MeV are listed in Tabled) These 
numbers are the average of the three measurements. The variation in the 7-ray 
flux between the three detectors with a distance of 50 cm from each other is 
about 3%. It is worthwhile mentioning that the contribution to the 7-ray flux 
from radon is estimated to be less than 1% for a radon level as high as 500 
Bq/m 3 , which is an example of the radon level measured at Homestake |23j . 
Therefore, the difference in the radon levels can not explain the difference seen 
in the measured 7-ray fluxes. The difference in the 7-ray fluxes measured for 
different levels corresponds to the variation of the radioactivity levels in the 
rocks near the experimental sites. The variation in the radioactivity levels is 
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Table 1 

The integrated 7-ray flux at the surface administration building, underground on 
the 800- ft, 2000- ft and 4550-ft levels. The quoted errors are the largest differences 
between the average of the three detectors and the individuals. 







The measured 


7-ray flux (cm 2 s 1 ) 








E>0.1 MeV 


E>1 MeV 


E>2 MeV 


E>3 MeV 




Surface 


1.56±0.05 


(4.63±0.14)xl0 -1 


(5.52±0.17)xl0~ 2 


(1.09±0.03)xl0- 


-3 


800 ft 


2.65±0.08 


(7.97±0.24)xl0 _1 


(9.49±0.28)xl0~ 2 


(4.81±0.20)xl(T 


-4 


2000 ft 


3.42±0.10 


(1.04±0.03) 


(1.26±0.04)xl0 _1 


(7.05±0.21)xl0~ 


-4 


4550 ft 


2.16±0.06 


(6.32±0.19)xl0- 1 


(9.64±0.29)xl0" 2 


(6.01±0.18)xl0- 


-4 



associated to the geologic formations for different levels as demonstrated via 
the Monte Carlo simulation described in Section 3. 



3 Simulated 7-ray flux for different levels at the Homestake Mine 

The measured 7 rays are primarily from the decay of actinide elements in 
the surrounding rocks. The contamination by the actinide elements and the 
rock composition must be taken into account in the Monte Carlo simulation. 
Since the highest 7-ray energy from natural radioactivity is 2.615 MeV and 
the attenuation length of such a 7 ray in the Homestake rock is about 9.4 cm, 
a rock thickness of 300 cm is sufficient to simulate the 7 rays from the rock. 
On the other hand, the attenuation length of a 7 ray with 2.615 MeV energy 
in air is about 21 meters. Therefore, we define the simulation geometry as 
below (Figure [2]): 1) the sensitive detection areas are two experimental halls 
10 mx 10 mx 3 m in dimension; 2) these two halls are separated by a 100 mx 
3 mx 3 m drift in the middle; 3) both the hall and the drift are filled with 
standard air; and 4) the outer space of the experimental hall and the drift are 
surrounded by a 3 m thickness of rock. 

A Homestake rock sample 278-2 [23] is applied in our simulation with the 
chemical composition given in Table [2j The chemical composition varies from 
location to location at different levels. Twenty rock samples near the surface 
were collected and analyzed [23] • The rock chemical composition affects the 
density and therefore the attenuation of the 7 rays through the rock. The 
largest variation in chemical composition seen in the samples was used to 
estimate the variation in the 7-ray fluxes. This variation results in a maximum 
uncertainty of about 10% in the simulated 7-ray fluxes. It is clear that the 
chemical composition cannot account for the variation seen in the 7-ray fluxes 
at different levels. 



5 




Fig. 2. Sketch of the geometry in the simulation showing two experimental halls 
intersected by a drift surrounded by 3 m of rock. 

Table 2 

The chemical composition of the Homestake rock sample used in the simulation |24] . 



Rock Component 
(Sample 278-2) 


Composition 
(% weight) 


Si0 2 


43.7±0.4 


TiC-2 


1.22±0.01 


A1 2 3 


13.6±0.1 


FeO 


12.7±0.1 


MnO 


0.13±0.01 


MgO 


7.0±0.1 


CaO 


7.9±0.1 


Na 2 


2.87±0.03 


K 2 


0.21±0.002 


P 2 5 


0.07±0.001 


H 2 


10.7±0.01 



The major contributions to the measured 7 rays are from U, Th, K and 
their daughters. The 7-ray decay ratios and energies from the decay chain of 
these isotopes were studied in detail in Ref . [28] . The simulation is performed 
by using the GEANT4 simulation package [29]. The 7 rays in a few MeV 
range can be simulated in an accuracy of less than 1% [2H] . The dominated 
uncertainty of about 10% comes from the variation of the chemical composition 
in the rock. The input 7 rays are uniformly generated in the rock with the 
emitted energies from the decay chains. We summarize the simulated 7-ray 
fluxes in Table |3j 

The simulated 7-ray flux is quoted in the number of 7 rays per ppm per cm 2 
per second. The ppm stands for parts per million describing the concentration 
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Table 3. The 7-ray flux in the simulated experimental hall from natural radioactivity in the rock. The rows indicate the energy of 7 ray 
in MeV size bins. The columns represent the further subdivision of the data in the MeV bins into 0.1 MeV bins. The uncertainty in the 
simulation is about 10%. 



Ranj 


s;e (MeV ) 


n 1 
U — .1 


1 
.1 — .z 


.1 — .O 


.3 - .4 


.4 -.5 


.5 - .6 


.0 — . 1 


7 Q 
. 1 — .O 


c 

.0 — .y 


n 1 

.y — 1 


Sum 


Hi™ 


oource 








7-ray flux (ppm 1 cm 2 s 














238 tt 
u 


q i P 9 


1 4p 1 


U.tJC Z, 


4.1e-2 


1.6e-2 


l.le-2 




7 fip ^ 


^ rip ^ 


^ 9p ^ 


A Dp 1 





232 Th 


3.5e-2 


5.5e-2 


3.3e-2 


1.2e-2 


6.7e-3 


1.2e-2 


3.0e-3 


5.3e-3 


2.7e-3 


7.5e-3 


1.7e-l 




40 K 


6.6e-2 


1.4e-l 


7.1e-2 


3.4e-2 


2.3e-2 


1.9e-2 


1.5e-2 


1.2e-2 


1.2e-2 


l.le-2 


4.0e-l 




238 TJ 


5.7e-3 


l.le-2 


7.4e-3 


6.1e-3 


4.1e-3 


1.7e-3 


1.0e-2 


1.4e-2 


1.5e-3 


2.2e-4 


6.2e-2 


1 


232 Th 


5.0e-3 


1.4e-3 


5.8e-4 


5.1e-4 


4.3e-4 


1.4e-3 


8.0e-4 


2.2e-4 


2.9e-4 


1.4e-4 


l.le-2 




40 K 


1.2e-2 


1.3e-2 


9.4e-3 


8.7e-3 


1.6e-l 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


2.0e-l 




238 TJ 


0.0e+00 


1.5e-3 


3.5e-3 


0.0e+00 


4.4e-4 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


5.5e-3 


> 2 


232 Th 


2.9e-4 


6.5e-4 


2.2e-4 


4.3e-4 


1.4e-4 


7.2e-4 


9.3e-3 


0.0e+00 


0.0e+00 


0.0e+00 


1.2e-2 




40 R 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 



Table 4 

The three main rock types and the contamination of 238 U, 232 Th and 40 K in the 
rocks [25] as measured from a number of samples. 







238TJ 


232 Th 


nat-K 




Average 


3.43 ppm 


7.11 ppm 


2.66 % 




Standard Deviation 

L> \J CJLA.L vlCXl LI J J V 1 < 1 U l^yll 


3 73 nnm 

tX . 1 U XX IJ XXX 




1.05 % 


Poorman Formation 


Number of samples 


7 


7 


141 




Maximum Value 


11.7 ppm 


18 ppm 


4.98 % 




lVTl 1"! 1 TYI 1 1 TT1 Vfl 1 1 1 P 


D 1 nDTTl 


ppm 


0.04 % 




Average 


1.51 ppm 


7.38 ppm 


0.96 % 




Standard Deviation 


ddtti 


4 04 nnm 

-±- * Lypiii 


0.91 % 


Homestake Formation 


Number of samples 


24 


24 


263 




Maximum Value 


3.5 ppm 


15 ppm 


4.31 % 




Minimum Value 


4 nnm 

w » [J k/111 


1 9 nnm 

X . i_y lj [Jill 


8 3e-3 % 




Average 


7.45 ppm 


30.5 ppm 


3.32 % 




Standard Deviation 


3.65 ppm 


16.5 ppm 


0.95 % 


Ellison Formation 


Number of samples 


2 


2 


55 




Maximum Value 


11.1 ppm 


47 ppm 


6.22 % 




Minimum Value 


3.8 ppm 


14 ppm 


1.16 % 



of radioisotopes in the rock, which allows us to calculate the 7-ray flux with 
the known concentrations of 238 U, 232 Th and 40 K. There are three main rock 
formations at the Homestake Mine: Poorman Formation, Homestake Forma- 
tion, and Ellison Formation [21~f25] . Table H] shows these three types of rock. 
Taking the average values of the concentrations, we obtain the 7-ray flux $ 
(E > 0.1 MeV) = 3.9±0.4, 2.3±0.2, and 9.5±1.0 cm^s" 1 for the Poorman, 
Homestake, and Ellision Formations, respectively. In addition, rhyolite intru- 
sive rocks occur as thin dikes in the Homestake rock [21] at different levels. 
The 238 U, 232 Th, and 40 K contents in the rhyolite dikes are higher than these 
in the above three formations. The counted rhyolite samples show a maximum 
of 9.4 ppm of 238 U, 12.2 ppm of 232 Th, and 3.98 % of K [21]. However, the 
rhyolite dikes do not make up large volumes in the underground [2*1] . 

To compare with the measured 7-ray fluxes at the 800-ft, 2000-ft, and 4550-ft 
levels, one must know the combination of the rock formation for these levels. A 
detailed geological study is under way. Using the preliminary knowledge from 
Refs. [21|2rj] , we assume that the rock near the measurement locations at the 
800-ft level is mainly dominated by the Poorman (80%) and Homestake (20%) 
Formations. This combination predicts a 7-ray flux (E>0.1 MeV) of 2.6±0.78 
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cm _2 s _1 , which is in good agreement with the measured 2.65±0.80 cm~ 2 s~ 1 at 
this level. The rock close to the measurements at the 2000-ft level constitute 
mainly Homestake (80%) and Ellison (20%) Formations. The predicted 7-ray 
flux (E>0.1 MeV) for the 2000-ft level is about 3.7±1.1 cnrV 1 . Again, it 
is in a good agreement with the measured flux of 3.42±1.0 cm _2 s _1 . The 
rock adjacent to the measurements at the 4550- ft level is in the Homestake 
Formation. The predicted 7-ray flux (E>0.1 MeV) of 2.3±0.7 cm _2 s _1 agrees 
with the measured 2.16±0.7 cm _2 s _1 very well. The variation of the 7-ray 
flux at different locations on the same level depends on the variation of the 
rock formations and the radioactivity levels in the rock. This variation can 
be as large as 30% as seen in the measurements. This is to say that the 7- 
ray flux must be measured in the experimental area where a low-background 
experiment is to be located. 



4 The simulated neutron fluxes for different levels 

Deep underground neutrons come from three primary sources: 1) the (a, n) 
neutrons induced by the a-decays in 238 U and 232 Th decay chains in the rock, 
2) the spontaneous fission decay of 238 U, and 3) the muon-induced neutrons. 
The muon-induced neutrons have a strong dependence on the depth [27]. In 
general the fluxes of the cosmogenic neutrons are one order, two orders, and 
three orders of magnitude smaller than that of the radiogenic neutrons at the 
800-ft level, the 2000-ft level, and the 4550-ft level, respectively. However, the 
cosmogenic neutron energy spectrum is much harder [27] and they are beyond 
the scope of this work. Neutrons generate 7 rays through inelastic scattering 
and capture on nuclei. These 7 rays have a higher energy than the 7 rays from 
radioactive decay described in the above sections. We simulate the neutron 
propagation in the rock with the same geometry that is used to simulate the 
7 rays described in Fig. [2J The mean free path of the radiogenic neutrons 
is about 10 cm in the rock. These neutrons can be substantially attenuated 
in the 3 m rock. Therefore, the contribution of the neutrons beyond 3 m is 
negligible. 

4-1 Neutrons from (a,n) Reactions 

The decays of 238 U, 232 Th, and their daughters produce a particles, 7 rays, 
and /? particles. The a particles with energies on the order of a few MeV range 
are capable of generating neutrons through (a,n) reactions in the rocks. The 
thickness of the rock in the simulation is 3 meters. 

The thick target yield of (a,n) neutrons for an a-particle projectile is calcu- 
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lated using the equations in Refs. [8f9] as below: 



r R N» r Eo (x-AE) 

Y t = l nME)^^l ^dE, (2) 

where R refers to a range of the a-particle in the materials. is the number of 
atoms per unit volume of the z-th element, Oi is the microscopic [a, n) reaction 
cross-section for an a-particle with the energy E, and Eq stands for the initial 
energy, S™ is the mass stopping power of the i-th element, A, is the atomic 
mass of the z-th element, and Na is Avogadro's number. 

Under the assumption of secular equilibrium, the 232 Th decay chain yields six 
a particles and the 238 U decay chain produces eight a particles with various 
energies E. The net neutron yield from the decay chains of 232 Th and 238 U 
can be determined by the sum of the individual yields induced by each a, 
weighted by the branching ratio for each element and weighted by the mass 
ratio in the rock. The energy attenuation of a-particles in the medium is the 
dominant process for the thick target hypothesis. Under the assumption that 
the incident flux of a-particles with energy of Ej is invariant until the energy 
is attenuated to zero, for target element i, the differential spectra of neutron 
yield can be expressed as: 



_ Na R*(E 3 ) da(E Q ,E n ) ^ 

~ jU STiEj) Jo dE a dba ' {6} 

where iVj is the total atomic number of the rock and <& a (Ej) is the flux of an a- 
particle with the specific energy Ej. R a (Ej) refers to the «-particle production 
rate for the channel with the energy Ej from an actinide element decay. If we 
consider the specific activity of an actinide with the concentration per ppm 
per gram per year, then 

R a ( Ei ) = 

Where A a stands for the mass number of the actinide element, specifically 
238 U or 232 Th here, Bj represents the a-particle branching ratio for the certain 
energy decay channel Ej, and tu2 is the half life of the decay. 

The thick target yield for a compound is discussed in Refs. [8f9] . In order 
to simplify the calculation, the compound is assumed to be a homogeneous 
mixture and Bragg's law of additivity for stopping power cross sections holds 
for the compound. Based on these assumptions and Eq.Q, the yield for a 
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compound Y c can be written as 



Y c = y £w i Y i 



(4) 



where Wi is the mass weight of element i in the compound. 

The cross section in Eq.©, e.g., J Ej M ^f n) dE a , is obtained from TALYS 
nuclear reaction code [30] in which the cross sections of neutron production 
for all possible reaction channels are calculated. The flux of a-particles is ob- 
tained by combining the production rate of a particles with the corresponding 
mass stopping power in the target. The mass stopping power for specific en- 
ergies is calculated by our simulation described in Ref. [31] and the ASTAR 
program [32]. We use the decay chains of selected isotopes in Ref. [33] for a- 
particle emission from 238 U and 232 Th decays. Only decays with visible energies 
larger than 0.1 MeV or branching ratio more than 0.5% are included. 

The overall uncertainty in the calculation of the net neutron yield is justified 
to be less than 20%. The sources of the uncertainty are the (a,n) reaction cross 
sections and the calculation of a stopping power in rock. The justification is 
done by comparing the net neutron yield from this work to similar calculations 
done for SNO [H] and Gran Sasso [TTj . 

4-2 Neutrons from Spontaneous Fission 

Naturally, neutrons can also be produced by spontaneous fission of the ra- 
dioactive elements, e.g., 238 U, 235 U, and 232 Th. The spontaneous fission of 238 U 
produces the most significant number of fission neutrons due to its shorter 
fission half-lives. The other fission neutron rates are two orders of magnitude 
smaller. Therefore, we consider the fission neutrons primarily from 238 U. The 
spectrum of the fission neutrons follows the Watt spectrum. 



where we adopt the Watt spectrum parameter a = 0.7124 MeV and b = 5.6405 
MeV -1 in the calculation. The production rate of spontaneous fission of 238 U 
is 0.52±0.02 neutrons ppm~ 1 g~ 1 year~ 1 in the rock [IT]. The uncertainty is 
from the measured 238 U decay constant for spontaneous fission [12J. 

4-3 The Simulated Neutron Flux and Neutron-Induced ^-Ray Flux 

The neutron flux and neutron-induced 7-ray flux in the experimental hall is 
predicted by the GEANT4 simulation package. The (a,n) induced neutrons 



N(E) = Cexp(-£/a) sinh(&£) 



1/2 



(5) 
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are isotropically generated in the rock with the energy spectrum given above 
for two different sources of actinide elements. Eq.Q is used to generate neu- 
trons produced by 238 U fission decay. 

Table and Table M show the simulated differential energy spectra of neu- 
trons induced by 238 U and 232 Th decay chains. These two Tables provide con- 
ventional values to be used in the simulation for predicting neutron-induced 
background for individual experiments if one knows the radioactivity levels in 
the rock that surrounds the detector. 

Tables [7J and [H] summarize the integrated neutron flux from (a,n) reactions 
and fission decays and the 7-ray flux induced by neutrons. The neutron flux 
is reported in terms of thermal neutrons, slow neutrons and fast neutrons. 
Note that the neutron fluxes reported in Table [7J are for the neutrons with 
kinetic energy greater than 0.1 MeV before they transport and scatter into 
the experimental hall. This means the fluxes of the thermal and slow neutrons 
and neutron-induced 7 rays are underestimated. 

Using the results from Tables [7J and [HJ the neutron and neutron induced 7- 
ray fluxes can be estimated by multiplying the flux with the radioactivity of 
actinide elements in the rock. There are core samples counted by Ref. [20] 
for different levels at the Homestake Mine. Utilizing the highest concentration 
of radioactivity counted for the 4850-ft level, i.e., 238 U: 0.55 ppm, 232 Th: 0.3 
ppm and 40 K: 2.21%; the prediction for the 4850-ft level is shown in Table |9j 
For the 7400-ft level, the highest measured concentrations [20] are 238 U: 0.49 
ppm, 232 Th: 0.20 ppm and K: 0.57%; giving a neutron yield and 7-ray flux in 
Table ED 



5 Discussion and conclusion 



Using three identical Nal detectors, we have investigated the 7-ray fluxes at 
the Homestake Mine at several levels that are currently accessible. A simple 
algorithm is developed to convert the measured 7-ray rates into a flux. The 
measured 7-ray fluxes are compared to a Monte Carlo simulation by applying 
known knowledge about the rock chemical composition, rock formation, and 
radioactivity levels. A good agreement is achieved between the measurements 
and Monte Carlo simulation. The Monte Carlo simulation is then applied to 
predict the neutron fluxes and 7-ray fluxes induced by natural radioactivity 
for different levels. These results are important for the initial set of physics 
experiments to be performed at the Sanford Lab. The predictions are useful 
to the planned low-background experiments at DUSEL. 
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Table 5. Neutron flux in the simulated experimental hall via (a, n) reaction due to natural radioactivity using the Homestake rock 
chemical composition. The rows indicate the energy of neutron in MeV size bins. The columns represent the further subdivision of the 



data in the MeV bins into 0.1 MeV bins. The uncertainty in 


the calculation is about 


20% that 


is justified 


in the text. 






Ranj; 


p (MeV) 


- .1 


.1 - .2 


.2 - .3 


.3 - .4 


.4 -.5 


.5 - .6 


fi 7 


7 8 
. ( — .o 


8 Q 

.o — .y 


Q 1 

.y — i 


Sum 


E n 


Source 








Neutron Yield (ppm 1 cm ' 


s ) 













238 JJ 


1.8e-6 


7.5e-8 


5.3e-8 


3.9e-8 


3.0e-8 


3.3e-8 


o.ue-o 


Q /Ip 8 


9 /lo 8 


1 7p 8 

1. 1 e-o 


9 9p p. 


232 Th 


6.7e-7 


2.6e-8 


2.0e-8 


1.4e-8 


9.9e-9 


1.2e-8 


i.oe-o 


1 ^p 8 

i .oe-o 


i ric 8 
i.ue-o 


fi Ao Q 

u.^e-y 


7 Qp 7 

i .ye- i 


1 


238 TJ 


1.3e-8 


1.9e-8 


1.8e-8 


1.7e-8 


1.7e-8 


1.9e-8 


1 fip 8 

i.oe-o 


1 % 8 

i.ze-o 


1 1 p 8 

i. re-o 


8 7p Q 

o. i e-y 


1 ^p 7 

r.oe- i 


232 Th 


5.0e-9 


7.9e-9 


7.2e-9 


7.4e-9 


5.6e-9 


6.5e-9 


t; 9 P q 
o.ze-y 


^ Dp Q 

o.ue-y 


7 fip Q 


q 7p q 
o. t e-y 


fi 1p 8 
u. ±e-o 




238 TJ 


1.3e-8 


7.7e-9 


l.le-8 


1.2e-8 


1.4e-8 


9.8e-9 


Q fip Q 

y.oe-y 


8 1 p Q 

o. ie-y 


fi 9p q 
o.ze-y 


fi 1 p Q 

u. ie-y 


Q 7p 8 

y. i e-o 


2 


232 Th 


4.9e-9 


4.1e-9 


6.9e-9 


4.3e-9 


4.9e-9 


4.9e-9 


4.ue-y 


o.ye-y 


o.ze-y 


9 Ac Q 
z.4e-y 


A Ap 8 




238 JJ 


6.8e-9 


7.5e-9 


6.7e-9 


3.0e-9 


3.1e-9 


1.8e-9 


1 fip Q 

i.oe-y 


o.oe-y 


z. ie-y 


9 7p q 

z. 1 e-y 


1 8p 8 

o.oe-o 


3 


232 Th 


2.7e-9 


2.9e-9 


2.0e-9 


1.2e-9 


1.8e-9 


1.6e-9 


i. ie-y 


1 1 p Q 

i . re-y 


1 ^P Q 

i.oe-y 


1 7p Q 
i . i e-y 


i 7p c 
i. i e-o 




238 JJ 


2.1e-9 


1.8e-9 


1.6e-9 


1.2e-9 


1.3e-9 


1.8e-9 


o.ye-iu 


1 fip_Q 

i.oe-y 


o.ye-iu 


1 9 D Q 

i.ze-y 


1 /Ip 8 


1 


232 Th 


1.0e-9 


1.6e-9 


l.le-9 


1.2e-9 


1.9e-9 


9.6e-10 


fi 4p 1 n 


c; Q P in 
o.ye-iu 


t; q p in 
o.ye-iu 


A in 


1 Dp 8 

i.ue-o 




238 JJ 


1.0e-9 


1.0e-9 


8.9e-10 


5.9e-10 


3.0e-10 


4.4e-10 


i .oe- ru 


8 Qp 1 n 
o.ye-iu 


tc Qp in 
o.ye-iu 


n noj_nn 
u.ue-t-uu 


o.ye-y 


5 


232 Th 


7.0e-10 


5.9e-10 


3.8e-10 


4.8e-10 


2.1e-10 


3.8e-10 


fi zip in 
u.^e-iu 


1 Q p in 
^.oe- ru 


1 Rp in 
i.ue-iu 


9 1 p 1 n 
z. ie- iu 


A 9p Q 


6 


238 JJ 


8.9e-10 


5.9e-10 


3.0e-10 


3.0e-10 


1.5e-10 


4.4e-10 


i.oe-iu 


n n c i nn 
u.ue-t-uu 


i.oe-iu 


1 Cnin 

i.oe-iu 


Q 1 p Q 

o. ie-y 


232^ 


l.le-10 


5.4e-10 


3.8e-10 


2.1e-10 


2.1e-10 


1.6e-10 


i . i e- 1 u 






Zi . i e- 1 u 


1 Qp Q 

i . ye-y 


7 


238 JJ 


1.5e-10 


3.0e-10 


1.5e-10 


0.0e+00 


4.4e-10 


0.0e+00 


A Ap 1 n 


u.ue~i-uu 


u.ue-t-uu 


q rip in 
o.ue-iu 


1 So Q 

i.oe-y 


232 Th 


l.le-10 


0.0e+00 


0.0e+00 


l.le-10 


0.0e+00 


0.0e+00 


l.le-10 


l.le-10 


5.4e-ll 


5.4e-ll 


5.4e-10 




238 jj 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


> 8 


232 Th 


5.4e-ll 


l.le-10 


5.4e-ll 


0.0e+00 


0.0e+00 


1.6e-10 


5.4e-ll 


0.0e+00 


5.4e-ll 


5.4e-ll 


5.4e-10 



Table 6. Neutron flux (ppm _1 cm~ 2 s~ 1 ) in the simulated experimental hall via 238 U fission decay due to natural radioactivity using 
the Homestake chemical composition. The rows indicate the energy of neutron in MeV size bins. The columns represent the further 
subdivision of the data in the MeV bins into 0.1 MeV bins. The uncertainty of about 3% that is stated in the text. 



E n (MeV) 


- .1 


.1 - .2 


.2 - .3 


.3 - .4 


.4- .5 


.5 - .6 


.6 - .7 


.7- .8 


.8 - .9 


.9 - 1 


Sum 





9.0e-7 


3.4e-8 


2.4e-8 


1.7e-8 


1.4e-8 


1.8e-8 


1.7e-8 


1.7e-8 


1.3e-8 


8.0e-9 


l.le-6 


1 


8.2e-9 


1.0e-8 


7.7e-9 


6.4e-9 


1.0e-8 


6.8e-9 


6.7e-9 


7.9e-9 


5.5e-9 


4.7e-9 


7.4e-8 


2 


5.1e-9 


4.7e-9 


7.7e-9 


4.7e-9 


4.5e-9 


3.7e-9 


3.8e-9 


3.5e-9 


2.9e-9 


2.4e-9 


4.3e-8 


3 


1.9e-9 


2.9e-9 


2.2e-9 


1.3e-9 


1.3e-9 


7.0e-10 


9.7e-10 


1.3e-9 


8.4e-10 


1.3e-9 


1.5e-8 


4 


1.2e-9 


9.0e-10 


l.le-9 


3.5e-10 


9.7e-10 


7.7e-10 


3.5e-10 


6.3e-10 


6.3e-10 


4.2e-10 


7.3e-9 


5 


3.5e-10 


3.5e-10 


3.5e-10 


2.8e-10 


2.8e-10 


5.6e-10 


1.4e-10 


4.9e-10 


0.0e+00 


0.0e+00 


2.8e-9 


6 


7.0e-ll 


0.0e+00 


7.0e-ll 


0.0e+00 


0.0e+00 


0.0e+00 


2.1e-10 


0.0e+00 


7.0e-ll 


1.4e-10 


5.6e-10 


7 


7.0e-ll 


1.4e-10 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


2.1e-10 


> 8 


0.0e+00 


7.0e-ll 


7.0e-ll 


7.0e-ll 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


0.0e+00 


2.1e-10 



Table 7 

The neutron flux induced by 238 U and 232 Th radioactivity in the simulated experi- 
mental hall. The quoted errors are justified in the text. 





Thermal Neutron 


Slow Neutron 


Fast Neutron 


Source 


E n <1 eV 


E n in [1 eV, 0.1 MeV] 


En >0.1 MeV 




(ppm cm _2 s ) 


(ppm cm _2 s ) 


(ppm cm _2 s —1 ) 


238 U (a,n) 


(1.17±0.23)xl0~ 6 


(6.79±1.36)xl(T 7 


(6.46±1.29)xl0~ 7 


232 Th (a,n) 


(4.19±0.84)xl0~ 7 


(2.48±0.50)xl0~ 7 


(2.64±0.53)xl0~ 7 


238 U fission 


(5.73±0.16)xl0~ 7 


(3.26±0.09)xl0~ 7 


(3.04±0.09)xl0~ 7 



Table 8 

The 7-ray flux induced by neutrons in the simulated experimental hall. The quoted 
errors are justified in the text. 

Source 7-ray flux (ppm cm~ 2 s ) 

238 U (a,n) (8.47±1.69)xl(T 7 

fission (3.92±0.11)xl(r 7 

232 Th (a,n) (3.03±0.61) x 1(T 7 

Table 9 

Predicted 7-ray and neutron fluxes for the 4850- ft and 7400- ft levels. The quoted 
errors are justified in the text. 



Depth 7-ray flux (cm 2 s 1 ) 
238 U and 232 Th 40 K 
4850 ft 0.32±0.10 1.46±0.44 
7400 ft 0.27±0.08 0.38±0.11 



Neutron flux ( cm 2 s 1 ) 

(2.3±0.8)xl0~ 6 
(2.0±0.7)xl0~ 6 
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